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Abstract. 

The radium atom is a promising system for studying parity and time invariance 
violating weak interactions. However, available experimental spectroscopic data for 
radium is insufficient for designing an optimal experimental setup. We calculate the 
energy levels and transition amplitudes for radium states of significant interest. Forty 
states corresponding to all possible configurations consisting of the 7s, 7p and 6d single- 
electron states as well as the states of the 7s8s, 7s8p and 7s7d configurations have been 
calculated. The energies of ten of these states corresponding to the 6d 2 , 7s8s, 7p 2 , and 
Qd7p configurations are not known from experiment. Calculations for barium are used 
to control the accuracy. 



PACS numbers: 31.25.Eb,31.25.Jf,32.70.Cs 
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1. Introduction 

Studying parity (P) and time (T) invariance violating effects in atoms is a way of 
searching for new physics beyond the standard model (see, e.g. [1]). These effects 
are strongly enhanced in radium atom due to high value of the nuclear charge Z, and 
specific features of the nuclear and electron structure j2JIHllIj- The atomic electric dipole 
moment induced by the T,P-violating nuclear forces and P-violating effects produced 
by the nuclear anapole moment are enhanced 3 orders of magnitude in comparison with 
previous experiments (the detailed comparison and complete list of references may be 
found e.g. in the review pQ). Preparations for the measurements are currently in 
progress at Argonne j3] and Groningen [SJ Ej • 

Detailed knowledge of the positions of the lowest states of an atom as well as 
transition probabilities between them is important for the design of cooling and trapping 
schemes and for estimation of the enhancement of the P and T-odd effects. Energy 
spectrum of Ra was first measured by Rasmussen [B] in 1934. Interpretation of his 
data was corrected by Russell [§] also in 1934. Compilation by Moore ^U] based on 
these two works contains about forty energy levels of radium. There were few more 
experimental works on radiums studying Rydberg states [TT], hyperfine structure and 
isotope shift [T2*l ITB*] . nuclear magnetic moments [T3], etc. In the most recent work 
by Sielzo et al the lifetime and position of the 3 P° state of Ra has been measured. 
The result for the energy is in excellent agreement with early data by Russell [9J and 
Rasmussen [B]. 

There were some doubts inspired by theoretical work of Bieroh et a/^H] among 
experimentalists working with radium on whether the data presented by Rasmussen jH] 
and Russell [Oj were reliable and accurate. This disagreement between theory 
and experiment motivated our previous calculations ^7]. The calculations strongly 
favored experimental data. However, the strongest evidence of the correctness of the 
experimental data came from recent success in trapping of radium atoms at Argonne. 
Corresponding paper which would include among other things new data on experimental 
resolution to the D-state location is to be submitted soon [THj . 

Excellent agreement between theory and experiment for radium and its lighter 
analog barium [T7] allows us to address next problem - gaps in experimental data for 
radium. In particular, it is important to know the positions of the states corresponding 
to the 6d 2 configuration. The locations of these levels are important when considering 
the possibility of laser-cooling and trapping Ba or Ra in the metastable 6s5d 3 D 3 or 
7s6d 3 D 3 state, respectively. This would be a useful alternative to the relatively slow 
and leaky transitions available from the ground 1 S'o state. In particular, the 6s5d 3 D% - 
5dQp 3 F 4 transition in barium and the 7s6d 3 _D 3 - Qd7p 3 F 4 transition in radium could 
provide a fast and closed cycling transition. However, the data for barium [THj indicates 
that the energies of the 5c? 2 configuration lie very low, between the Qs6p and the 5d6p 
configuration, and therefore provide an undesirable leak channel. It is reasonable to 
expect that the energies of the 6d 2 configuration of radium also lie pretty low. This would 
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limit the cooling and trapping schemes for radium causing leaking of some transitions 
into the states of the 6d 2 configuration [18 . 

The main task of our previous paper was to prove that the experimental data 
was correct. Therefore we calculated only energy levels known from experiment. In 
present paper we extend the calculations to include all states of the lowest configurations 
of radium. We calculate energy levels and lifetimes of forty states of the 7s 2 , 7s7p, 7s6d, 
7s8s, 6d7p, 6d 2 , 7s8p, 7p 2 and 7s7d configurations. For 19 lowest states we also present 
detailed data on electric dipole transition amplitudes. Similar calculations for barium 
are used to control the accuracy of the calculations. 



2. Method of calculations and results for barium 

The method of calculations has been described in detail in our previous works |20| l2T] 
E21 Ulj • Here we repeat its main points most relevant to present calculations. 

The calculations are done in the V N ~ 2 approximation [21 which means that initial 
Hartree-Fock procedure is done for a double ionized ion, with two valence electrons 
removed. This approach has many advantages. It simplifies the inclusion of the core- 
valence correlations by avoiding the so called subtraction diagrams[2IIJ EI]- This in 
turn allows one to go beyond second-order of the many-body perturbation theory 
(MBPT) in treating core-valence correlations. Inclusion of the higher-order core- valence 
correlations significantly improves the accuracy of the results [2UE2]- Also, using V N ~ 2 
approximation makes calculations for a positive ion and for a neutral atom very similar 
providing more opportunities for the control of the accuracy. One more advantage is 
that atomic core is independent on the state of valence electrons. Ground and excited 
states are treated equally which is important for calculating energy intervals. 

Single-electron Hamiltonian for a valence electron has the form 

h^ho + ti, (1) 

where ho is the relativistic Hartree-Fock Hamiltonian: 

Ze 2 

h = cap + (3- l)mc 2 - — + V N ~ 2 , (2) 

r 

and Ei is the correlation potential operator which represents correlation interaction of 
a valence electron with the core. Calculations for a positive ion are done by solving the 
equation 

(hi - e v )i) v = 0, (3) 

where e v and ip v are the energy and wave function of a valence electron. Both e v and ip v 
include the effect of core-valence correlations and the wave functions ip v are often called 
Brueckner orbitals to distinguish them from Hartree-Fock orbitals which do not include 
correlations. 
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The effective Hamiltonian for a neutral two-electron atom is the sum of two 
single-electron Hamiltonians plus an operator representing interaction between valence 
electrons: 

H eS = h 1 (r 1 ) + h 1 (r 2 ) + h 2 (r 1 ,r 2 ). (4) 

Interaction between valence electrons is the sum of Coulomb interaction and correlation 
correction operator £2: 

e 2 

h 2 = 1 1 + E 2 (n,r 2 ), (5) 

l r l - r 2| 

£2 represents screening of Coulomb interaction between valence electrons by core 
electrons. 

We use standard configuration interaction (CI) technique to solve the Schrodinger 
equation for two-electron valence states. Two-electron wave function for the valence 
electrons \l/ has a form of expansion over single-determinant wave functions 

^ = J2ci®i(n,r 2 ). (6) 

i 

$i are constructed from the single-electron valence basis states calculated in the V N ~ 2 
potential 

$i(ri, r 2 ) = -^(^ a (n)^(r 2 ) - ^b(n)ij a (r 2 )). (7) 

Coefficients q as well as two-electron energies are found by solving matrix eigenvalue 
problem 

{Htf - E)X = 0, (8) 

where Hff = (^{H^l^) and X = { Cl , c 2 , . . . , c n }. 

The most complicated part of the calculations is calculation of the correlation 
correction operators Si and S 2 . We use MBPT and Feynman diagram technique to 
do the calculations. MBPT expansion for £ starts from second order. Inclusion of the 
second order operators Y*^ and £^ effective Hamiltonian (0J) accounts for most of 
the core-valence correlations. However, further improvement is still possible if higher- 
order correlations are included into Si. We do this the same way as for a single valence 
electron atoms . Two dominating classes of the higher-order diagrams are included 
into £1 by applying Feynman diagram technique to the part of £1 which corresponds to 
direct Coulomb interaction. These two classes are (a) screening of Coulomb interaction 
between valence and core electrons by other core electrons, and (b) interaction between 
an electron excited from the core and the hole in the core created by this excitation [23J. 
The effect of screening of Coulomb interaction in exchange diagrams is imitated by 
introducing screening factors fk into each Coulomb line. We assume that screening 
factors fk depend only on the multipolarity of the Coulomb interaction k. It turns out 
that the values of fk vary very little from atom to atom and the same values can be 
used for all atoms of the first and second columns of the periodic table: 

fo = 0.72, /1 = 0.62, f 2 = 0.83, f 3 = 0.89, / 4 = 0.94, / 5 = 1.0, .... 
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Calculations show that for atoms like Ba and Ra accurate treatment of Ei is more 
important than that of S 2 . Therefore we calculate £2 i n second order of MBPT only. 

One needs a complete set of single-electron states to calculate E and for construction 
of two-electron basis states (JJJ) for the CI calculations. We use the same basis in both 
cases. It is constructed using 5-spline technique [SHEET We use 50 S-splines of order 
7 in a cavity of radius R max = 40a#, where a& is Bohr radius. Single-electron basis 
orbitals in each partial wave are constructed as linear combination of 50 5-splines 

50 

Mr) = Y,b?B t (r). (9) 

i=l 

Coefficients bf are found from the condition that ip a is an eigenstate of the Hartree-Fork 
Hamiltonian ho (J2J)- 

The effect of inclusion of second and higher-order £ into effective Hamiltonian for 
two-electron valence states of Ba and Ra was studied in detail in our previous paper |17j . 
It was also suggested there that the best results can be obtained if fitting parameters 
are introduced before for each partial wave. The values of these parameters for 
Ba found from fitting experimental energies of the 6s 2 , Qs6p and 6s5d configurations 
are X s = 1.0032,A P = 1.0046 and A^ = 0.9164. Note that we keep the same fitting 
parameters for £ Pl/2 and £ P3/2 as well as for £d 3/2 and £d 5/2 - We do this to avoid 
false contribution to the fine structure. Fitting of the energies imitate the effects of 
higher-order correlations, incompleteness of the basis set, Breit and QED corrections. 

Final results for Ba are presented in Tabled The results for twelve states of the 6s 2 , 
6s5d, 6s6p and 5d6p configurations are the same as in our previous work. However we 
present now 27 more states, including states of very important 5<i 2 configuration. Note 
that corresponding energies absent in Moore book JU| and we use recent compilation by 
Curry ^H] instead. Parameter A in the Table is the difference between experimental and 
theoretical energies (A = E expt — E ca i c ). The agreement between theory and experiment 
is extremely good in most of cases. The largest difference is for the 5c? 2 1 S state. It 
is 723 cm" 1 or 2.7%. Note however that experimental value for this state came from a 
different source than all other data and has the largest uncertainty (see (TSJ for details). 
There is a chance that the experimental value is incorrect. The only other large difference 
is for the 5g! 2 state. It is 409 cm -1 or 1.8%. For other states of the 5g? 2 configuration 
the difference between theory and experiment is about 1% or smaller. For most of other 
states the difference is just small fraction of a per cent. 

Energy levels of barium where calculated by many authors before 123 I2H1 12H 
EDI CHj - The scope of the present work does not allow us to cite all these results. 
Comprehensive review of previous calculations for Ba is a big task while our present 
consideration serves very specific and limited purpose. We just want to demonstrate 
that our method work very well for Ba, therefore we can expect the results of similar 
quality for Ra which has similar electron structure. 

In Table d we also present the values of observed and calculated g-factors. Non- 
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Table 1. Energies and lifetimes of lower states of barium 
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Table 2. Experimental and theoretical transition probabilities for barium 
Lower Upper Transition probability (s -1 ) 
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relativistic (NR) values are given by 

J(J + 1)-L(L + 1) + S(S + 1) 

9NR = 1 + 2J(JT1) ' (10) 

where J is total momentum of the atom, L is angular momentum and S is spin. 
Comparing calculated values of ^-factors with observed and non-relativistic values is 
useful for identification of the states. 

In Table ^ we also present calculated values of lifetimes of all considered states. 
Only electric dipole (El) transitions were included in the calculations. Therefore, we 
don't present lifetimes of the long living states which can only decay via magnetic dipole 
(Ml) or electric quadrupole (E2) transitions. 

We calculate El transition amplitudes between states \l/ a and ^ using the 
expression 

A(El) ab = Y,c?c b J ($ l \\d z + 5V N - 2 \\® 3 ), (11) 

where d = — er is the electric dipole operator, 5V N ~ 2 is the correction to the self- 
consistent potential of the atomic core due to the electric field of the photon. The 
term with 5V N ~ 2 accounts for the so called RPA (random-phase approximation) or core 
polarization correction. The functions are two-electron basis states (J7J) and c, are 
expansion coefficients for states \l/ a and over basis states <3>; as in (jOJ. 

Expression (jllj) is approximate. It includes dominating contributions to the El 
amplitudes but doesn't take into account some small corrections. A detailed discussion 
of different contributions into matrix elements between many-electron wave functions 
can be found e.g. in Ref. [32J. In terms of that paper expression ([TTj) corresponds to 
the leading contribution to the effective amplitude (Arpa, see Eq. (22) of Ref. [32J). It 
accounts for configuration interaction, core-valence correlations and core polarization 
effects. Next, the so called subtraction contribution (Asbt) does not exist in present 
calculations since we use the V N ~ 2 approximation. Subtraction terms appear only if 
Hartree-Fock procedure includes valence electrons. They account for the difference 
between Hartree-Fock potential and potential of the core in the CI Hamiltonian. In the 
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V N ~ 2 approximation for a two valence electrons atom these two potentials are identical. 
The terms not included into (jllj) are: the two-particle correction (Arp),the self-energy 
correction (A a ), structure radiation and normalization corrections (see j3Z| for details). 

Using expression (fTTj) gives satisfactory accuracy for most of the cases. However 
neglecting other contributions for small amplitudes may lead to some instability of the 
results. This is especially true for small amplitudes which vanish in the non-relativistic 
limit (AS* > 0, AL > 1). Present calculations give only rough estimation of the values 
of these amplitudes. However, it doesn't have much effect on lifetimes since lifetimes 
dominate by strong transitions with large amplitudes. 

Typical accuracy of the calculations for strong transitions is illustrated by the data 
in Table El Here we compare some calculated transition probabilities for Ba with the 
most accurate experimental data. The probability of the El transition from state i to 
a lower state j is (atomic units) 

3. Results for radium 

The results of calculations for energies, ^-factors and lifetimes of forty lowest states of 
radium are presented in Table El Energies are compared with available experimental 
data. Calculations follow the same procedure as for barium. The only difference is in 
values of rescaling parameters for correlation potential S. Fitting of the experimental 
energies of the 7s 2 , 7s7p and 7s6d configurations leads to the following values of the 
rescaling parameters: X s = 1.0021,A P = 1.0053 and = 0.9327. These values are very 
close to similar values for barium (see above). The Coulomb integrals and correlation 
corrections in electronic analogues (e.g. Ba and Ra) usually have approximately the 
same values. This fact may be used to extract unaccounted higher correlation corrections 
from Ba and improve our predictions for unknown energy levels in Ra. Indeed, the 
differences between theory and experiment for similar states of radium and barium 
are very close at least for lower states. This is in spite of different order of levels, 
about 2.5 times difference in fine structure intervals (spin-orbit interaction increases 
~ Z 2 ) and some difference in fitting parameters for the correlation potential. This 
means that the difference between theory and experiment for barium can be used to 
improve the predicted positions of those states of radium for which experimental data 
is absent. Column Extrap. in Table El presents energies of radium corrected using the 
difference between theory and experiment for barium. States where experimental data is 
available illustrate that the procedure leads to systematic improvement of the agreement 
between theory and experiment for lower states of radium. For states where there is no 
experimental data extrapolated values give better prediction of the energies than just 
ab initio calculations. 

Note that this procedure doesn't work for higher states. This is because saturation 
of the basis in the CI calculations rapidly deteriorates with the increase of the excitation 
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Table 3. Energies and lifetimes of lower states of radium 
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30982 




30573 




1.00 


1.05 


150 ns 


7s8p 


3po 





31085.88 


31008 


78 


31135 


49 


0.00 


0.00 


76 ns 




3po 


1 


31563.29 


30695 


868 






1.50 


1.07 


20 ns 




3po 


2 


31874.44 


31778 


96 


31909 


35 


1.50 


1.44 


57 ns 


7p 2 


3p 







29840 












21 ns 




3p 


1 


31248.61 


31365 


-116 






1.50 


1.49 


26 ns 




3p 


2 


32941.13 


33180 


-239 






1.50 


1.21 


42 ns 


7s7d 


3 D 


1 


32000.82 


31895 


106 


31968 


33 


0.50 


0.51 


18 ns 




3 D 


2 


31993.41 


31902 


91 


31981 


12 


1.17 


1.16 


19 ns 




3 D 


3 


32197.28 


32068 


129 


32155 


42 


1.33 


1.33 


21 ns 


7p 2 


X D 


2 


32214.84 


32205 


10 






1.00 


1.20 


29 ns 


6d7p 


3 D o 


1 


32229.97 


32090 


140 






0.50 


0.84 


21 ns 




3 D o 


2 


32506.59 


32436 


71 






1.17 


1.17 


13 ns 




3 D o 


3 


33197.46 


33169 


28 






1.33 


1.17 


21 ns 


7s8p 


lpo 


1 


32857.68 


31446 


1412 






1.00 


1.16 


34 ns 


6d7p 


3po 





33782.41 


33809 


-27 






0.00 


0.00 


10 ns 




3po 


1 


33823.70 


33837 


-13 






1.50 


1.40 


10 ns 




3po 


2 


34382.91 


34421 


-38 






1.50 


1.42 


11 ns 


6d 2 


x s 







33961 








0.00 


0.00 


150 ns 


6d7p 


lpo 


3 




34332 








1.00 


1.14 


25 ns 


6d 2 


x s 







35408 








0.00 


0.00 


30 ns 


6d7p 


lpo 


1 




36043 








1.00 


1.03 


38 ns 
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Table 4. El-transition amplitudes for 19 lowest states of radium 



Even state Odd state Amplitude Even state Odd state Amplitude 



7s 2 


% 


7s7p 


3po 


1.218 


7s6d 


X D 2 


6d7p 


1 D 2 


5.704 






7s7p 


lpo 


5.504 






Qd7p 




0.774 


7s8s 




7s7p 


3po 
r 1 


0.057 


Qd 2 


3 F 2 


7s7p 


3po 


0.542 






7s7p 


lpo 
r 1 


4.176 






7s7p 


lpo 


0.442 


7s6d 




7s7p 


3po 
r 


2.995 






7s7p 


3po 
r 2 


0.266 






7s7p 


3po 
r 1 


2.574 






6d7p 


3po 
r 2 


4.644 






7s7p 


lpo 


0.437 






6d7p 


u 2 


1.208 






7s7p 


3po 
r 2 


0.688 






6d7p 


3po 
r 3 


1.786 






6d7p 


3po 
r 2 


3.729 


Qd 2 


X D 2 


7s7p 


3po 


1.274 






6d7p 


If)" 
u 2 


1.394 






7s7p 


lpo 


1.023 


7s8s 




7s7p 


3po 

r o 


2.214 






7s7p 


3po 
r 2 


1.535 






7s7p 


3po 


3.890 






6d7p 




0.290 






7s7p 


lpo 


1.476 






6d7p 




3.259 






7s7p 


3po 
r 2 


6.075 






6d7p 




0.595 






6d7p 


3po 
r 2 


0.266 


7s6d 


3 D 3 


7s7p 


3po 
r 2 


6.340 






6d7p 


l v° 2 


3.584 






6d7p 


3po 
r 2 


0.107 


7s6d 


3 D 2 


7s7p 


3po 


4.382 






6d7p 




2.911 






7s7p 


lpo 
r 1 


0.813 






6d7p 


3po 
r 3 


3.064 






7s7p 


3po 
r 2 


2.605 






6d7p 


3po 
r 4 


5.885 






6d7p 


3po 
r 2 


2.946 


Qd 2 


3 F 3 


7s7p 


3po 
r 2 


0.190 






6d7p 


u 2 


0.168 






6d7p 




0.702 






6d7p 


3po 
r 3 


4.568 






6d7p 


1 D 2 


0.566 


7s6d 




7s7p 


3po 
r 1 


0.344 






Qd7p 


3po 
r 3 


5.672 






7s7p 


lpo 


3.189 






6d7p 


3 F^ 


1.597 






7s7p 


3po 
r 2 


0.510 


Qd 2 


3 F 4 


6d7p 




0.037 






6d7p 


3 F° 


2.856 






Qd7p 


3 F^ 


6.343 



energy. Since the energies of similar excited configurations of Ba and Ra are significantly 
different the effect of incompleteness of the basis is different too. 

Experimental data for g-factors of radium is not available. However, comparing 
calculated and non-relativistic values of ^-factors indicates that the L — S scheme still 
works very well for the most of the lower states of Ra and can be unambiguously used to 
name the states. The L — S scheme breaks higher in the spectrum due to the combination 
of relativistic effects and configuration mixing. For example, as can be seen from Table El 
states 7s8p 3 P?, Qd7p 3 D? and 7s8p 1 P° are strongly mixed. The calculated ^-factors of 
each of these states deviate significantly from the non-relativistic values. This makes it 
difficult to identify the states. Also, strong configuration mixing is probably the reason 
for poor agreement between theory and experiment for the energies of these states. 



Calculation of energy levels and transition amplitudes for barium and radium 



11 



Lifetimes of the states presented in Table|3]were calculated using Eqs. (jTTjl and (fT2|) 
for all possible electric dipole transitions from a given state to lower states. This involves 
270 El-transition amplitudes. It is impractical to present all of them in a table. However, 
for considering different trapping and cooling schemes it is important to know transition 
probabilities between different pairs of states rather than just lifetimes. Therefore we 
present in Table |U 52 amplitudes between 19 lowest states of radium. This data should 
be sufficient in most of cases. More data is available from authors on request. Note that 
the values of small amplitudes which vanish in non-relativistic limit (AS > 0, AL > 1) 
should be considered as rough estimation only (see discussion in previous section). 
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